Absence of local thermal equilibrium in two models of heat conduction 
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A crucial assumption in the conventional description of thermal conduction is the existence of local thermal equilibrium. We 
test this assumption in two simple models of heat conduction. Our first model is a linear chain of planar spins with nearest 
neighbour couplings, and the second model is that of a Lorentz gas. We look at the steady state of the system when the two 
ends are connected to heat baths at temperatures Tj and T2. If 21 = T2, the system reaches thermal equilibrium. If Ti ^ T2, 
there is a heat current through the system, but there is no local thermal equilibrium. This is true even in the limit of large 
system size, when the heat current goes to zero. We argue that this is due to the existence of an infinity of local conservation 
laws in their dynamics. 



PACS numbers: 05.70.Ln, 05.20.-y, 44.10. +i, 75.10.Hk 



Heat conduction provides a rather elementary example 
of nonequilibrium steady states (NSS). Here the driving 
force on the system is an externally applied temperature 
gradient. Our present understanding of this in terms of 
nonequilibrium statistical mechanics is not entirely sat- 
isfactory. In the traditional treatment of heat conduc- 
tion, a very important assumption is of a rather quick 
( in atomic time scales) establishment of local thermo- 
dynamical equilibrium (LTE) [Q. This allows one to 
define thermodynamical quantities such as temperature, 
pressure etc. locally. One then writes evolution equa- 
tions for the slow change in time in these quantities in 
terms of their small spatial gradients. Thus, in heat con- 
duction, LTE implies that we can define a temperature 
field T(x), that varies slowly in space. In linear response 
theory, the local heat current density J{x) is given by 
J = -K{T)VT (Fourier's law). The well-known Kubo 
formula determines the thermal conductivity coefficient 
K(T) in terms of the time-dependent equilibrium corre- 
lations of the system j2j . 

However, the question of the necessary (or suffi- 
cient) conditions for the fast establishment or non- 
establishment of LTE has not been investigated much 
so far. This is what we study in this paper, in two simple 
models. The first is a one-dimensional spin model evolv- 
ing with a Markovian stochastic dynamics. The second 
model is that of a d-dimcnsional (d > 1) Lorentz gas of 
nonintcracting particles scattering elastically with a set 
of randomly placed obstacles {e.g spheres) of finite den- 
sity. In both these models, we find the unexpected result 
that the NSS does not show LTE, even in the limit of 
very small thermal gradients. Thus, one cannot define a 
local temperature, and the conventional linear response 
theory does not work. We argue that a sufficient condi- 
tion for LTE not to occur is the existence of an infinity 
of locally conserved quantities in the dynamics. 



It is obvious that existence of extra conservation laws 
in an isolated classical mechanical system implies break- 
down of ergodicity, and hence failure to reach thermody- 
namical equilibrium. Our results are nontrivial in that 
we show that this occurs even in the presence of coupling 
to heat baths which leads to breakdown of the conserva- 
tion laws at the boundaries of the system. Also, the role 
of conservation laws in stochastically evolving system is 
somewhat different from that in deterministic evolution 
||, and not so well-studied. 

Several microscopic models of heat conduction have 
been studied in the past. The simple model of a harmonic 
lattice, evolving with classical mechanical equations of 
motion, and coupled to two heat baths at different tem- 
peratures Ti and T2 at opposite ends of the lattice, is 
known to show anomalous conduction It is found 
that the heat current J across the sample remains finite 
as the size of the system L goes to infinity, for a fixed 
Ti and T 2 . Simulations of harmonic chains with disorder 
show that J decreases to zero for large L, but only as 
L^ 1 ! 2 Several models with nonlinear couplings have 
been studied numerically: the FPU chain with quartic 
potential Q , the Toda lattice JjJ , the so-called ding-a-ling 
model ||, the Frenkel-Kontorova model |9| etc. Most of 
these models are one-dimensional, and all do not show 
the expected J ~ (7\ — T 2 )/L behavior. The uniqueness 
of the steady state has been proved only for some spe- 
cial models of heat baths, and only if the temperature 
difference is sufficiently small pOj . In this respect mod- 
els with stochastic dynamics have been more successful. 
These typically work with local energy-conserving moves, 
and often involve introduction of additional degrees of 
freedom. Creutz has used an algorithm with Maxwell 
demons to simulate heat conduction in the Ising model 
|llf . Lebowitz and Spohn studied heat conduction in the 
Lorentz model, and showed that in the Boltzmann-Grad 
limit of large number of scatterers of very small size, one 
recovers the Fourier law ( J ~ {T x -T 2 )/L) @. 
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We start with a precise definition of our first model. 
We consider a linear chain of L planar spins. The spin at 
site i (1 < i < L) of the lattice is specified by the angle 0j, 

< 9 < 2ir. The spins interact with nearest neighbors 
by ferromagnetic coupling K. The Hamiltonian of the 
system is given by 

H = ~K <x*(Pi-0i) (1) 

<i,j> 

where the sum is over all nearest neighbors. 

The dynamics is the following: Suppose the instanta- 
neous local field at a site % due to coupling of the neigh- 
bors is in the direction fa, and the spin at the site is 

01 = <pi + 5&i. Then the transition 50i — > — S9i does not 
change the energy of the system. We assume that such 
spin flips occur at all sites stochastically with a constant 
rate (which may be chosen to be 1). The boundary spins 
i = 1,L are connected to heat baths. For these spins, 
the flip rates are the following: A boundary spin Oi can 
change to any value 9\ with a rate a if the energy change 
AE < 0, and with a rate ae~ AE/T , if AE > 0; where T 
is the temperature of the heat bath. 

Note that in this model energy is conserved exactly 
away from the boundaries of the lattice. In the absence 
of any coupling to heat baths, this dynamics has been 
studied by one of us earlier |13fj . It was found that, 
in general, the dynamics is not ergodic, and the phase 
space breaks up into disconnected sectors. We now show 
that the coupling to heat bath allows transitions between 
these sectors, and makes the system fully ergodic. 

Since the rate matrix in our problem satisfies the de- 
tailed balance condition, it is sufficient to show that all 
configurations can be reached from any initial configura- 
tion. Consider first the case when a single spin on the 
lattice, say 0q, is coupled to a heat bath at temperature 
T. Then, 9q can be made to take any value, since it 
can exchange energy with the bath, keeping other spins 
unchanged. Next consider any spin, 9\ that is nearest 
neighbour to Qq. Then it is easy to see that flipping 0\, 
then changing 9q by a small amount SOq, then flipping 
9\ again, then changing 0$ back again to its initial value 
keeps all spins unchanged, except for 9\ which changes 
by an amount proportional to &9§. 

We can use a similar argument to change a neighbour 
of Oi, and so on. Thus, there exists an allowed sequence 
of spin flips by which we can rotate any single spin by 
an infinitesimal amount. By making many such rotations 
we can reach any spin configuration. This completes the 
proof of ergodicity. Note that the proof is valid in all 
dimensions d and also if more than one spin is in contact 
with the same heat bath ( same T). 

In one dimension the equilibrium properties of the XY 
model are easily determined. We define A9i = 6>; + i — 
Then clearly the A9i are independent random variables. 
For a mesoscopic description, we define coarse-grained 



densities by averaging the corresponding microscopic 
quantity over length scales I, with 1 << I « L. Let 
u^ n \x) be the coarse-grained density corresponding to 
cos n (A6'i) for i lying in a neighborhood of the point x. 
In equilibrium, at inverse temperature (5, there is no de- 
pendence on x, and one gets 

(n) _ J cos n (A9)e l3Kcos ^d{A9) 

U i - J e/ 3ifco S (A0) d ( A6 /) ( 2 ) 

Eliminating [3K from these equations we can express 
vS n \ for n > 2, as explicit nonlinear functions of vP-', 
i.e. = FW(«W). 

We consider now the case when the two ends are con- 
nected to different heat baths, at temperatures T\ and 
T 2 . We show that, in the NSS, LTE is not achieved. 

In this dynamics, whenever a spin flip occurs, the val- 
ues of A9 on the two adjacent bonds get interchanged. 
This means that the densities are locally conserved 
for all positive integers n. For each of these quantities 
one can write a corresponding current. The current from 
the (x — l)th bond to the xth bond is given by 

jW(a?) = (- < cos"(A^_i) > + < cos"(A6y >) 

= -Vu^(x) (3) 

In the steady state, in one dimension, J^ n >{x) must be 
independent of x, and hence vary linearly across the 
chain for all n. Their value at the two ends are deter- 
mined by the temperatures at the ends. Then, eliminat- 
ing the coordinate x, u^ n \x) is expressible as a linear 
function of u^(x), for all n. But, as shown above, in 
equilibrium, , for n > 2 are nonlinear functions of 
u^K Thus, we have constructed a simple model in which 
heat current in the steady state is locally proportional 
to the gradient of energy density, but there is no LTE 
(Fig. 1). This result is true even in the limit of system 
size L — > oo, when the heat current through the lattice 
becomes infinitesimal. 

For a simple extension of the above result, consider a 
linear chain in which there is one bond with a different 
bond-strength K' . Then the quantities (n > 2) are 
not conserved at this 'bad' bond. In this case, the cur- 
rents take different values j[ and to the left and 
right of the bond respectively (for n > 2). In this case, 
there are no strict conservation laws, but are still 
linear functions of i/ 1 ) in each (left or right) half of the 
chain, and LTE is still not obtained. This conclusion has 
also been checked in simulations (Fig. 2). 

In higher dimensions, there are no known locally con- 
served quantities, other than energy, in the spin model. 
Hence, we expect LTE to be attained. Indeed, a numer- 
ical simulation of a 40 x 40 lattice with T\ = 1AK, and 
T-2 = OAK agrees fully with this conclusion. In this case, 
we checked that the observed value of < u' 2 ' (x) > at any 
point x in the NSS agrees very well with the value of 
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in the homogeneous equilibrium state having energy den- 
sity u^(x). Similarly, if we consider a dimerized chain 
with alternating coupling constants K and K' , then there 
are no conserved quantities other than the energy den- 
sity. In this case, the results of our simulations ( also 
shown in Fig. 1) confirm that the NSS does show LTE 
even in one dimension if K ^ K' . 

Our second example of absence of LTE in heat con- 
duction is provided by the model of the d-dimensional 
Lorentz gas studied earlier by Lebowitz and Spohn Jl2| . 
The model describes a gas of non-interacting point parti- 
cles moving in a in box, undergoing elastic scattering by 
a random assembly of fixed obstacles of arbitrary shape 
(Fig. 3). In this model, all collisions with the obsta- 
cles are elastic, and energy is conserved. But collisions 
between the particles and the two walls (at x = and 
x = L) are inelastic, and lead to the energy after colli- 
sion being thermalized corresponding to the temperature 
of the wall. 

Let p(E, x)dE be the average density of particles hav- 
ing energy between E and E + dE in a small vol- 
ume centered at point x. Then all moments p^ n '[x) = 
J E n p(E, x)dE of this distribution function are locally 
conserved. 

Now consider the effect of coupling it to two differ- 
ent reservoirs at temperatures T± and T^. Let p(x) be 
the probability that a randomly chosen particle in the 
small volume chosen near x was introduced at the left 
end. Then as collisions with scatterers do not change the 
energies, clearly p(E, x) is given by 

p{E, x) = P (x)p(E, x = 0) + (1 - p(x))p(E, x = L) (4) 

The x-dependence of p(E, x) comes only from the spa- 
tial dependence of p(x). As the linear combination of 
two maxwellians is not a maxwellian, and the distribu- 
tion must be a maxwellian in thermal equilibrium, it fol- 
lows that there is no thermal equilibrium in the Lorentz 
model jL4|. Also we expect the heat flow to be diffusive 
( J ~ in the limit L >> £, where I is the mean 

free path of the Lorentz particles. Note that an equation 
similar to Eq. (4) can be written down for the first model 
also. 

If we allow inelastic scattering with local energy con- 
servation ( each scatterer can store a small amount of 
energy, which can be exchanged with the scattered parti- 
cles ) , the infinity of conservation laws goes away, and we 
expect that then the model would show LTE, and normal 
conduction. 

To summarize, in this Letter, we studied two simple 
models of heat conduction, and showed that local thermal 
equilibrium is not reached as both the models have an 
infinity of locally conserved quantities. These counter- 
examples should help understand better the mechanism 
of local thermal equilibriation in nonequilibrium systems, 
in general. 
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FIG. 1. Plot of u^\x) versus u^ x \x) in the steady state 
of the linear chain (K = 1). The results for the dimerized 
chain (K' — 2) agree with the equilibrium curve. 




FIG. 3. A schematic representation of the Lorentz model. 
Particles move ballistically, and are scattered elastically, but 
in a random direction on collision with an obstacle. At the 
left and right boundaries, they are reflected and are given a 
new energy randomly from a distribution corresponding to 
two different temperatures Ti and T2 
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FIG. 2. Plot of versus when the middle bond 
has a different strength K' — 0.5 showing the different linear 
dependences in the left and right halves of the chain. For 
comparison, the case K' = 1, and the equilibrium curve are 
also shown. 
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